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1 Experiment, Observation and Astronomy 

Astronomers observe. 

Most other sciences permit experiment of some kind, where cause is ma- 
nipulated and effect monitored; as astronomers, however, we are denied this 
luxury and must learn what we can of the Universe by eavesdropping on Na- 
ture as she murmurs to herself. Denied the luxury of experiment, we rely on 
experience guided by intuition to interpret, and so bring order, to our observa- 
tions. Our experience is so far removed from the environment we are studying, 
however, that it is a poor guide to our imagination. 

An important lesson of astronomy is that the ways we see determine the 
things we know. We thought the Universe a different place when, before 
Galileo, we saw it only through the naked eye: the real failure of the Ptole- 
maic model was only manifest when we saw that the planets had disks and 
that Venus passed through a full set of phases. Each increase in sensitivity, 
spatial or temporal resolution, or spectral range has brought with it new puz- 
zles and surprises, and ultimately new insights and understanding. Progress in 
astronomy has historically been associated with new or improved observational 
perspectives. 

Thus, to compensate for the vast difference between our own and Natures 
laboratories, we have learned to rely on a multiplicity of observational perspec- 
tives: we view the Universe in as many different ways as we can in the hope 
that the multiplicity of distinct perspectives will lead us to a more reliable 
view of Nature's workings. 

When completed, the gravitational wave detectors now proposed or under 
construction will provide us with a perspective on the Universe fundamentally 
different from any we have come to know. With this new perspective comes the 
hope that new insights and understandings of Nature will emerge. The pro- 
posed acoustic detectors, with spherical geometries and operating at millikelvin 
temperatures, are particularly well-suited to study gravitational radiation in 
the 1-3 KHz band. In this brief report I review some sources of particular 
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interest for these detectors. 

2 The detector 

For the purpose of describing the sensitivity to astrophysical gravitational ra- 
diation snuxces of the next generation acoustic detectors, I assume an array 
of TIGA^dEI "spheres" operated together as a broadband detector. Each el- 
ement in the array has a different resonant frequency /o and a bandwdith 
A/o ~ Jo/10; over this bandwidth the one-sided noise power spectral density 
in either h+ or h x is constant. The resonant frequencies and bandwidths of 
the different array elements are arranged so that they adjoin with negligible 
overlap so that broadband detector composed of these separate elements covers 
the band from 850 Hz to 2650 Hz with a one-sided noise power spectral density 
of approximately 10~ 46 Hz _1 in either h+ or h x - Finally, outside of this band 
I assume that the noise power is effectively infinite. 

3 Lumpy neutron stars 

Among the more interesting gravitational radiation sources for a high frequency 
detector like .that described in section^ is a rapidly rotating, non-axisymmetric 
neutron starQ Setting aside pulsar spin-down over the course of an observation 
as well as the time-varying Doppler shift arising from Earth's rotation on its 
axis and revolution about the Sun, the radiation from such a lumpy neutron 
star, rotating about a principle moment of inertia with rotational period 1/ f m t, 
is 

4 *| _|_ cog^ i 2 

h+ = (27r/ rot ) e/ 3 cos 47r/ rot t (1) 

r 2 

4 2 

h x = — cosi (27r/ rot ) e/ 3 sin47r/ rot i (2) 
r 

where Ij are the neutron star's principal moments of inertia, rotation is about 
the I3 axis, I2 — I1 = and i is the inclination angle between the neutron star 
angular momentum and the detector line-of-sight. The mean square signal-to- 
noise (averaged over all i) in our hypothetical detector for an observation of 
duration T is 

In order to interpret the signal-to-noise ratio we can evaluate the false 
alarm rate: the probability that, in a fixed observation, noise alone will result 
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in a signal-to-noise ratio above a fixed threshold. A search for a continuous- 
wave source in a time series of duration T can achieve a frequency resolution no 
better than 1/T. The filtered detector output in each frequency "bin" of width 
A/ = 1/T is independent. The filtered output can be further decomposed 
into an in-phase and quadrature-phase component (relative to an arbitrary, 
monochromatic reference signal), which are also independent. The total signal- 
to-noise, which is the Pythagorean sum of the signal-to-noise of the in-phase 
and quadrature filtered components, is Rayleigh distributed: 

P(p) = pexp(-p 2 /2); (4) 

consequently, the probability that p < po is 

C(p ) - 1 - exp (-pl/2) (5) 

and the probability that, across all bins in a bandwidth /, all the signal-to- noise 
ratios are less than po is 

P( P < p ) = C(p ) fT . (6) 

The false alarm rate (i.e., the probability that at least one bin has p > po) is 
thus 1 — C(p)f T . Assuming (as we have in section ||) a detector bandwidth 
of 1150 Hz and an observation period of 1 yr, p > 7.5 corresponds to a false 
alarm rate (odds that one or more of the 3.6 x 10 10 "bins" are above threshold 
owing to noise) of 2.%, or 1 false alarm every 50 y -1 . 

The. likely maximum strength of the neutron star crust provides an upper 
bouncEra on e of 10 -6 . Beyond this upper bound, theory provides no guidance 
on e. Observation, on the other hand, does. Gravitational radiation emis- 
sion contributes to the spin-down of a rapidly rotating neutron star; thus, the 
observed spin-down rate of any millisecond pulsar limits e for that pulsar. Ob- 
served millisecond pulsars have very low spin-down rates and correspondingly 
low upper-bounds on e: no observed millisecond pulsar admits an e greater 
than a few times 10~ 8 (Blair, private communication). 

Were it known that for some part of the millisecond neutron star popu- 
lation e did take on a reasonably large value then one coul/Luse the expected 
ms neutron star space density (n ~ 10 -6 pc~ 3 in the disku'a) to estimate the 
expected number with signal-to-noise greater than some threshold pq. Let's do 
that. Suppose that 1% of ms neutron stars have e ~ 10~ 7 and rotational fre- 
quencies in the range 100 < / rot /Hz < 1000. Assume that these are distributed 
in frequency so that dn/df oc f~ 2 , where n is the number space density. Take 
the galactic disk radius to be 15 Kpc, with Sol 8 Kpc from the galactic center, 
and the ms pulsar disk scaleheight H to be 300 pc (the same as for LMXBs). 
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Restricting attention to those sources with p > 7.5 in a 1 y observation we can 
expect 

• 27 sources out to the near edge of the galactic disk in the band 850 Hz < 
/ gw < 1 KHz; 

• 240 sources in the band lKHz < / gw < 1.8 KHz; 

• 350 sources throughout the galaxy in the band 1.8 KHz < / gw < 2 KHz. 

A precessing neutron star will also radiate gravitationally in a manner very 
similar to that of a non-axisymmctric neutron star rotating about a principal 
axis. For the small wobble angles that might be expected in a realistic source, 
the radiation amplitudes are typically much sniauer than for the case studied 
here. The interested reader may find details inl§E£l. 



4 Coalescing compact binaries 

Coalescing binary neutron star systems are a significant gravitational wave 
source for the ground-based inteferometric gravitational wave detectors like 
LIGO and VIRGO. These systems radiate most of their energy at low frequen- 
cies, making them less suited as sources for the proposed acoustic detectors. 
Nevertheless, it is instructive to consider how observable compact binary sys- 
tems are to the proposed acoustic detectors. 

The amplitude signal-to-naise ratio for binary coalescence in an omni- 
directional detector is given byllnla 

'= 8 (£)(i^r ec(/ ™->' < ? > 

where ro is a characteristic distance that depends only on the source waveform 
and the detector noise power spectrum, oIl is the binary's luminosity distance, 
depends on the angular orientation of the binary with respect to the line 
of sight to the detector and is of order unity, ( describes the fraction of the 
detector bandwidth that is filled by radiation from the inspiraling binary sys- 
tem, and / CO ai is the orbital frequency (not the gravitational wave frequency) 
where the inspiral_ends and the coalescence begins. For binaries with equal 
mass components^ 

/ coal ~710H^T^ (8) 
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where M is the binary system's total mass, and for the detector described in 
section ^, 



f 4/3 ^ f < /min/2 ' 

C(/) = I '^-fS' 3 if /-»<// 2< /— ( 9 ) 

{ 1 if / > /max/2 



/min 


= 850 Hz, 




(10) 


/max 


= 2650 Hz, 


and 


(11) 


ro 


= 28Mpc ( 


'l0- 46 Hz\ 

s n )■ 


(12) 



For neutron star binaries (mi = m,2 = 1.4 Mq) £ is approximately 1/2. 

The overlap function £ taken together with the relationship between the 
system mass and the coalescence frequency / coa i shows immediately that the 
inspiral radiation from binaries with total mass greater than about 4.7 M Q 
will not leave any imprint on the proposed detectors (however, see the discus- 
sion of black hole formation in section |E|). Similarly, the value of ro for the 
proposed detector suggests that, in any event, binaries not much further than 
the Virgo cluster wilLJpe visible. Assuming that the rate density of neutron 
star binary inspiraolij is 3y -1 at 200 Mpc, neutron star binary coalescence 
within 20 Mpc is not expected more frequently than once every 330 y. The 
inspiral of "conventional" compact binaries is thus not a source of interest for 
the advanced proposed acoustic detectors. 

What of unconventional binary systems, however? Recent results from the 
MACHO collaboration suggest that perhaps 50% of the galactic halo mass is 
in dark objects with m ~ 0.3M Q [L 6 J Assuming that 

• a fraction w ~ 50% of the halo mass ( Mhaio greater than or approximately 
equal to 6 x 10 11 Mq) is in MACHOs, 

• all of the MACHOs are mass m ~ 0.3 M Q neutron stars or black holes, 

• a fraction x ~ 30% of the MACHOs are bound in symmetric binaries, 

• this binary population was formed at the same time that the galaxy was 
formed (age T ~ 10 10 y), and 

• the MACHO binary population is coalescing today at a steady rate, 
then the rate of MACHO binary coalescence in the halo is given by 



If we take the distance to the typical halo binary to be 50Kpc, then the 
amplitude signal-to-noise in our antenna array of a typical halo inspiral is 

p~510. (14) 

Thus, under these (optimistic!) assumptions, a TIGA antenna array would see 
virtually all halo binary inspirals: 15 strong (p ~ 500!) events per year. 

5 Black hole formation 

When a binary whose total mass M is greater than the maximum neutron 
star mass coalesces the likely outcome is a mass M black hole. The time 
dependence of the spacetime strain owing to the excitation of fundamental 
quadrupole mode of the resulting blackhole is, at late times, given by 

m(t) = e ^ ftlQ sin (2tt/£) (15) 

where Q and / are related to the black hole mass and angular momentum. 

In general there are five fundamental quadrupole modes as well as a succes- 
sion of "overtones." The fundamental modes are degenerate if the black hole 
is non-rotating (Schwarzschild) and are split otherwise. The overtones are at 
higher frequency and are also more strongly damped than the fundamental 
modes. To give a rough estimate of the detectability of the gravitational ra- 
diation from the late stages of black hole formation, let's ignore both rotation 
and the overtones and assume, forpCpnvenience, that black hole excitation is 
concentrated in a single mode; therO 

/ ~ 12. KHz (^), (16) 
Q ~ 2. (17) 

For Schwarzschild black holes, the proposed detector is sensitive to for- 
mation of black holes with mass between 4.5 and 14.1 M Q . For comparison, 
the mass of Cygnus X-l is estimated to be in the range 5-10 M Q , LMC X-3 is 
about 10 M Q , and A0620-00 is about 3.8 M©; so, the formation of astrophysical 
black holes in our antenna array's "mass range" is reasonable. 

Assuming a quadrupole radiation pattern for the black hole we can relate 
the "efficiency" of the black hole formation (ratio of the total energy radiated 
to the black hole mass) to the rms (averaged over all orientations of the black 
hole with respect to the detector line-of-sight) strain induced in each element 
of our antenna array: 

h lms {t) = 2^~e—m(t) (18) 
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Numerical modeling!^ of Schwarzschild black hole head-on collisions give effi- 
ciencies of 10~ 4 . For collisions of rotating black holes with significant orbital 
angular momentum the efficiencies might be significantly higher — perhaps 
1% or greater. 

The mean square amplitude signal-to-noise ratio of Schjwarzschild black 
hole formation in our hypothetical detector is thus given byLH 

-= l e \/20Mpc\ 2 / M \ 3 /10- 46 Hz _1 \ , , 

Even under the pessimistic assumption that there is no increase in efficiency 
for the coalescence of orbiting black holes our antenna array is sensitive to, 
we can expect to observe a typical 10 M© black hole formation event with a 
signal-to-noise of several out to the Virgo cluster. 

The rate of black hole formation is entirely uncertain; however, most as- 
trophysicists see no reason why the same mechanisms that make neutron stars 
do not also make black holes at approximately the same rate. By our present 
understanding of formation mechanisms, this is not a high rate even at the 
distance of the Virgo cluster; so, evidence of black hole formation at any mea- 
surable rate would require a significant change in our understanding of stellar 
evolution. 



6 Supernovae 

Theoretical models of stellar core collapse, and the corresponding gravitational 
wave luminosity, have a long and checkered history: estimates of the gravita- 
tional wave luminosity have at different times ranged over more than four 
orders of magnitude. It is not simply the luminosity that is unknown: the 
waveforms themselves are also uncertain, leading to a further difficulty in es- 
timating the detectability of this source. Nevertheless, it is still possible to 
evaluate what is required of stellar core collapse in order that it be observable 
in our hypothetical detector. 

Suppose that the waveform from supernovae is given by 

h+ = 2M® a /+ro(t) (20) 
r 

r 

a and (3 are constants, /+ and / x are functions of the relative orientation of the 
source with respect to the detector, and m(t) is some function of time which 
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we leave undetermined for now. The power radiated into each polarization 
mode is given by 

E + = a 2 (/ 2 )M||m| 2 (22) 

E x = /? 2 (/ 2 )M||m| 2 , (23) 

where <> signifies an average over the sphere. 

Now assume that equal power is radiated into the two polarization modes. 
Then we can write a and (3 in terms of a single parameter e as 

a 2 = I (2 a) 

2M (/ 2 )/^|mP 1 > 

p2 = 2M (/ 2 )/^H 2 - (25) 
In terms of e the power radiated into the + and x polarization states is thus 

eM |m| 2 , , 

S + = ^ = 27W (26) 

Now let us return to consider the time dependence of the waveform m(t). 
Note that 

' dt\m(t)\ 2 = I d/(2^/) 2 |m(/)| 2 , (27) 



where m is the Fourier transform of m. Assume that there is equal radiated 
power in equal bandwidths out to a frequency /o; then |m| 2 / 2 is constant for 
/ < fo and zero for / > /o- The mean-square signal-to- noise ratio (averaged 
over the source orientation relative to the detector) in the detector described 
in section || is then 

where 

/min = max[/ min ,min(/ ,/min)] (29) 
/max = max[/ min ,min(/o,/ max )] ■ (30) 

For the detector described in section ||, / m ; n is 850 Hz and / max is 2650 Hz. If 
we assume that / m ; n < fo ~ 1 KHz < / max , then 

/ 2\ 19n e ^ 100Kpc V l0- 46 Hz- 1 10 3 Hz 850Hz(l-/ min // o ) ^ oi ^ 
{P)C " 42 %0-H r J S n f / min 0.15 



Current calculations (which, bear in mind, are still unable to successfully 
describe the supernova explosion) suggest .£- in the range 10 _9 -10 _8 Mq, with 
peak power in the 200-300 Hz bandE3E3 : EJ Thus, without a very optimistic 
efficiency and a substantially wider signal bandwidth, we can't expect to be 
able to observe supernovae much beyond our own galaxy, and certainly not 
out to the Virgo cluster. 

7 Conclusions 

There are several interesting opportunities for an array of millikclvin acoustic 
gravitational antennae covering the 1-3 KHz frequency band. In particular, we 
can expect that the radiation from population of rapidly rotating neutron stars 
with oblateness on order 1/10 that allowed by the crustal breaking strength 
would be observable throughout the galaxy in a one year observation. Local 
supernovae, neutron star binary inspiral, or black hole formation, while in- 
frequent, would be a serendipitous radiation source from our own galaxy or, 
perhaps, as far as the Virgo cluster. Finally, a very speculative source — co- 
alescence of a population of compact MACHOs in the galactic halo — would 
be observable with large signal-to-noise ratio. 
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